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ABSTRACT 

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained 
as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary 
exhibits strong C m /14650 emission and is classified as an 06.5 Iafc supergiant, while the secondary is an 06 Iaf supergiant. The radial-velocity 
variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence 
for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M\ sin 3 i = 78 ± 8 Mq 
and M 2 sin 3 i = 66 ± 7 Mq. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide 
an excellent test for atmospheric and evolutionary models. 
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1. Introduction 

Massive binary stars provide vital insights to our understand 
ing of massive-star evolution. This is primarily due to the ac- 
curacy with which their masses can be determined: an essen 
tial ingredient for understanding a wide range of stellar proper 
ties. Because of the additional constraints that can be placed on 
their age and evolution, these stars provide information on ini- 
tial masses, chemical mixing and mass-loss (Moffat 2008] |De 
Mink et al.|2009| >. In a broader context, they can then act as cru- 
cial calibration points for models of both stellar atmospheres and 
evolution. 

From the catalogue of bright stars in the Magellanic Clouds 
by |Feast et~aL] ( |1960[ l, R139 has a V-band magnitude of ~1Z 
making it one of the brightest objects in the 30 Doradus nebulaH 
Walb orn & Blades| ( [TW7] l noted R 139 as potentially one of the 
most massive stars in 30 Doradus, urging more detailed study. 

Multi -epoch spectroscopy of R139 was obtained as part of 
the campaign by Moffat ( 1989 1. His mean radial velocity from 
observations in 1982 showed an offset of ~100kms _1 compared 
to the mean velocity from earlier data. He noted R139 as a 
single-lined binary, with a tentative period of 52.7 d adopted 
from a number of possible fits to the data. |Schnurr et al.] ([2008b ) 
presented spectroscopy of R139 from three observing seasons 
(spanning 2001 to 2003). While noting that the system displayed 
'slightly variable radial velocity', it was concluded that R139 
was single, citing the relatively large uncertainties in Moffat's 
past work for the conflicting scenarios. 



R139 has now been observed as part of the VLT-FLAMES 
Tarantula Survej|^](VFTS), an ESO Large Programme which has 
obtained multi-epoch spectroscopy of over 800 massive stars in 
30 Doradus. A full overview of the survey is given by Evans et 
al. (Paper I, 201 1), here we report on the discovery of R139 as a 
massive double-lined binary. 

pe] 2. Observations 

The observations of R 139 are summarised in Table [T] The pri- 
mary dataset of the VTFS has been obtained with the Giraffe 
spectrograph using the 'Medusa' fibre-fed mode of the FLAMES 
instrument on the Very Large Telescope (VLT). Details of the re- 
ductions and observational strategy can be found in Paper I, in 
which R139 is catalogued as object VTFS 527. 

After the initial detection of binarity made from the 
FLAMES data, follow-up observations have been obtained on 
the 6.5m Magellan Clay Telescope with the MagE instrument, 
with X-Shooter on the VLT and also with FEROS on the 
MPG/ESO 2.2m telescope at La Silla^] All these follow-up ob- 
servations provide coverage across the entire visible spectrum, 
with a typical signal-to-noise ratio of order 150 - although the 
FEROS data must be degraded to R ~ 9000 to achieve this. 

3. Results 

3.1. Binary identification and spectral classification 



* Hubble Fellow 

** European Space Agency 

1 Other a liases of R1 39 include: Brey 86 jBreysacher| |l98l| , 
Parker 95 2 jParker||1993|l, BA T99-107 < |Breysacher et al.||1999) ,"and 
Selman 2 ( |Selman et al.|1999) . 



Feast et aI7| ( |1960 > described the spectrum of R 139 as exhibiting 
'O-type absorption plus weak W emission' . |Walborn & Blades] 

2 Observations obtained at the European Southern Observatory Very 
Large Telescope in programme 182.D-0222. 

3 FEROS observations obtained as part of programme 086.D-0997 
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Table 1. Observational epochs for R139 detailing; the instru- 
ment used, the resolving power (R), the date and the time elapsed 
since the first epoch (A HJD). 
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1 1 1M 1 L1111C11L / iJCLllllg 




HID +7400000 


A HJD 
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1 


PT AMP'S /T T?09 


7 nnn 
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2 


FT AMFS /T R02 




54748 8278 


0.046 


3 


FLAMES /LR02 




54749.7233 


0.946 


4 


FLAMES /HR15N 


16,000 


54749.7726 


0.996 


5 


FLAMES /LR03 


8,500 


54755.6882 


6.911 


6 


FLAMES /HR15N 




54810.6697 


61.893 


7 


FLAMES /LR03 




54810.7374 


61.961 


8 


FLAMES /LR03 




54810.7808 


62.004 


9 


FLAMES /LR03 




54810.8357 


62.059 


10 


FLAMES /LR02 




54837.6332 


88.856 


11 


FLAMES /LR02 




54868.5569 


119.780 


12 


FLAMES /LR02 




55112.8129 


364.036 


13 


MagE / 1" slit 


4,500 


55527.8017 


779.025 


14 


X-Shooter/0.5" slit 


9,100^ 


55580.5477 


831.771 


15 


X- Shooter 




55584.7076 


835.931 


16 


X-Shooter 




55588.6922 


839.915 


17 


FEROS 


48,000 


55604.6909 


855.936 


18 


FEROS 




55605.5613 


856.739 



The FLAMES observations were comprised of back-to-back observations so 
only the mid-exposure HJD is given. 

The resolving power quoted for X-shooter is that of the UV-arm only, which 
overlaps the LR02-03 regions of the FLAMES data. 



( |1997| l argued that R139 has very strong Of emission features, 
classifying it as 07 Iafp, while Schnurr et al. (2008b) adopted 
the spectral type WN9h::a. This ambiguity has likely arisen be- 
cause the emission features are the superposition from two sim- 
ilar stars. When observed at lower resolution, this would have 
falsely suggested enhanced emission or broadening of the lines. 

The increased resolution and time-sampling of the new data 
have revealed that many of R139's prominent emission and ab- 
sorption features separate into two distinct and similar compo- 
nents. This is best illustrated by contrasting epochs 11 and 16 in 
Figure[T] which show the minimum and maximum observed sep- 
arations respectively. The epochs which display well-separated 
components have allowed a precise classification of both com- 
ponents. The system consists of a more massive and more lu- 
minous primary, which is an 06.5 Iafc supergiant, and a slightly 
less luminous 06 Iaf companion. These spectral types are deter- 
mined from visual inspection of the He n / He i absorption line 
ratios between ,14200 / /14026 and also between ,14542 / ,14471, 
the latter of which can be seen in Figure [T] 




Fig. 1. Normalised spectra of R139 at epochs of minimum 
(#11) and maximum (#16) observed separation. The He i ,14471 
line profile in epoch #16 suffers from a slight nebular over- 
subtraction but does not prevent clear identification of the com- 
ponents. 




Fig. 2. Normalised spectra showing the complex N m and C in 
emission region. Both of the N in emission lines and the He n 
,14686 emission are comprised of two un-equal components, the 
stronger of which is associated with the primary. The C in emis- 
sion lines show no separation, but exhibit radial velocity shifts 
in the same direction as the primary. 



The Ifc classification of the primary arises from emission of 
C m ,14647 and ,14650-4652. This relatively rare feature led to 
the recent introduction of the Ofc category in the morphological 
framework used to classify O-type spectra ( Walborn et al. 2010] >. 
Figure [2] clearly illustrates that the N in emission lines separate 
into two distinct components, whereas the two C in emission 
lines show no evidence of separation, but merely shift between 
epochs in the same sense as the primary. The Ofc feature has 
previously been associated with the 05 spectral type ( Walborn] 
|et al.|2010) ; this therefore, is an interesting example of the phe- 
nomenon in a later-type star, albeit in the LMC. 

Where possible, the relative shift of the components was 
identified through the C m emission and also He i ,14922 emis- 
sion, which are only present in the primary. In the LR02 ob- 
servations, the Si iv ,14116 emission line was used as the main 
diagnostic for the relative shift of the components. The Struve- 



Sahade effect (e.g. Linder et al. 



glected in this system given its re 



2007) can most likely be ne- 



atively long period. 



3.2. Radial velocity analysis and lower mass limits 

A global x 1 fitting approach has been used to determine the ra- 
dial velocity shifts of the different epochs. This technique fits 
double Gaussian profiles to a number of lines: He i ,14026, Si iv 
,141 16, He ii ,14200, S iv ,14485, S iv ,14505, He n A4542, He n A 
4686 and He i ,14922. The fitting is performed simultaneously on 
all the observations, which ensures that consistent profile shapes 
are used, including at conjunction. This approach improves the 
disentangling of the contribution from each star for the data sets 
with limited phase coverage, but ignores the possibility of line 
profile variations. The formal errors on the measurements for 
each component are a few kms -1 . 

The mass ratio of the system can be found from the ratios of 
the primary and secondary radial velocities, and is independent 
of any other assumptions about the orbit ( |Rauw et al.|2000| >. For 
R139 the mass ratio is found to be: Mi/M 2 = 1.20 + 0.05. 

Period searches based on Fourier analysis of the measured 
radial velocity shifts were performed using the methods of 
|Gosset et al. ( 200 1 ) > : the dominant signal indicated a period of 
153.9 days. The corresponding orbital solution, which has an 
rms uncertainty of 6.2 km s -1 , is shown in Figure [3] (based on 
methods of Sana et aT||2006a i. The corresponding lower mass 



limits for the stars are: M\ sin 3 i = 78 + 8 Mq and M2 sin 3 i - 
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Fig. 3. The best-fit orbital solution from the measured radial ve- 
locities of the components, indicating a 153.9 day orbit. The 
closed circles denote the velocities of the primary and the open 
circles the secondary. 

66 + 7 Mq. Other orbital parameters are listed in Table [2] These 
are sensitive to the behaviour of the system around periastron: 
further observations of this stage in the orbit would allow confir- 
mation of these results. 



3.3. The luminosity of R1 39 

To estimate the luminosity of the system, model atmospheres 
were calculated with CMFGEN (|Hillier & Miller|1998|>, adopt- 



ing abundances from |Aspfund et aL (2005 1 and scaling them ap- 
propriately for the LMC. These were used to constrain the ef- 
fective temperature (T e g) consistent with: 1) the absence of N iv 
/14058 emission; 2) the presence of N m /14640 emission; 3) the 
intensity of the He i AAAll absorption. This gives an estimated 
T e s for both components of 34 ± 2 kK. 

The luminosity was then determined by matching optical 
and infrared photometry from |Selman et al.] ( |1999| l and 2MASS 
( Skrutskie et al.||2006| >. For this purpose, the visual extinction 
Ay = R X E(B - V) was determined for each model, based 
on the relation R = 1.12 x E(V - K)/E{B - V) + 0.02 from 



Fitzpatrick (1999i. The resulting luminosity for the composite 
system is log(L/L G ) = 6.4 ± 0.1 (with R in the range 3.4-3.9). 



4. Photometric Variability 

R139 was identified as showing slight photometric variability by 
Feitzinger & Isserstedt (1983). They observed a 0.3 magnitude 
dimming in the V band over a 25 day period. However, this was 
from only three observations taken with a wide (18") aperture. 



Table 2. The parameters associated with the best-fit orbital solu- 
tion. The errors quoted are the formal errors on the best-fit from 
the Fourier analysis, and therefore may not be fully representa- 
tive of the uncertainty in the parameter values. 

Property Best-fit value 

Period, P 153.9 ± 0.1 days 

Eccentricity, e 0.46 ± 0.02 

Argument of periastron, cj 106.9 ±5.0 deg 

Date of <p = (HJD - 2450000), T 6035.9 ±1.3 days 

Maximum velocity of primary, Kl 107.8 ± 3.8kms~' 

Maximum velocity of secondary, K2 127.0 ± 4.5 km s -1 

Projected semi-major axis for primary, ai sin i 290.6 ± 10.8 Rq 

Projected semi-major axis for secondary, a2 sin i 342.3 ± 12.8 Rq 
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Fig. 4. Differential V-band residuals for R139 compared to the 
mean of the five check stars, phased to the 153.8 day orbit. Open 
circles denote the five epochs when the seeing was in excess of 
270. The dotted lines indicate the expected mid-point of any pos- 
sible eclipses: <p = 0.39 and 0.98 for the eclipses near apastron 
and periastron respectively. 



An active component of the VTFS is photometric follow-up 
with the 2 m Faulkes Telescope South, which has been used to 
monitor seven fields in the 30 Dor region. The default mode of 
the camera is 2 x 2 binning of the CCD pixels, giving an effective 
pixel-scale on the sky of (J.'278. 

We have 54 V-band epochs for the relevant field, spanning 
an 18 month period starting in January 2009. The Faulkes data 
are reduced automatically following observations, but are not 
calibrated photometrically. Given the crowding in this field, we 
used apphot in iraf to obtain instrumental magnitudes of R 139 
from aperture photometry. Five 'check' stars of similar bright- 
ness were selected from the frames for comparison: R133, R137, 
R138, Mk 1 1, and R146. From these, differential residuals (AV) 
were calculated for R139 compared to the mean magnitude of 
the check star for each epoch. The deviation was found to be con- 
sistent with that calculated between the check stars themselves, 
indicating that R139 shows no photometric variability. These re- 
sults are shown in Figure |4] where the observations have been 
phased to the 153.9 day orbit. 

If the inclination (/) of the system was 90°, the maximum 
duration of eclipses near apastron and periastron has been cal- 
culated to be 7.9 and 2.9 days respectively. From the sampling 
of our photometric data, it is unlikely that such events would 
have gone undetected, see Figure [4] However, if the inclination 
is lower (80° < i < 86°), there is no eclipse near apastron and 
the periastron eclipse is shorter. Consequently, an intensive pho- 
tometric observing campaign is required near to periastron to 
conclusively determine if there is any evidence for an eclipse. 



5. X-rays 



R139 was detected by Portegies Zwart et al.| (|2002) as an X- 
ray source in the 30 Dor field observed with the Advanced 
CCD Imaging Spectrometer on the Chandra X-Ray Observatory. 
Further analysis of the Chandra data was carried out by 



Townsley et al. ( 2006 1 and more recently by |Guerrero & Chu 
(2008J. These studies found R139 to have a relatively low X- 
ray luminosity compared with other W-R stars in the region. 
Guerrero & Chu ( 2008 ) also considered data from the Rontgen 
Satellite (ROSAT), but it did not detect R139 owing to its lower 
sensitivity. 
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Fig. 5. Hertzsprung-Russell diagram showing mass estimates de- 
rived from luminosity fits to evolutionary tracks (Brott et al. 
201 1; Friedrich et al. in prep.). The red circle indicates the best 
fit for the total luminosity of the R139 system, while the blue 
squares show the fit for the luminosity of the two components 
based on the mass ratio and assuming the objects are coeval. 
From these it is possible to infer the initial masses, see text for 
further details. 

The X-ray luminosity and the bolometric luminosity of mas- 
sive O stars are linked by the relationship Lx ~ 10~ 6 '• L\, i 
(Sana et al. 2006b). Therefore, with a luminosity of log(L/L Q ) = 
6.4 for the combined system, an X-ray luminosity of 1.2 x 
10 33 ergs~' would be expected. This is slightly lower than 
Guerrero's result of 2.7 x 10 33 erg s in the 0.5 - 7.0keV range 
- even considering the possible 25% error in the detected count 
rate. This slight excess emission might be associated with X-rays 
generated through the interaction of the system's stellar winds. 
There is no evidence, however, for phase-dependent line profile 
variations, which would have also suggested colliding winds. 



6. Discussion 

Previous observations: In order to compare our result with the 
earlier work of Schnurr et al. ( 2008b I, who found an insignificant 
radial velocity shift, our LR02 observations were degraded to the 
same resolving power as Schnurr's (R ~ 1000) and a number of 
lines were fit with a single Gaussian function. The radial velocity 
variation was found to be only 10.3 km s , while the FWHM of 
the profiles varied by around 40kms~ 1 . This suggests that, even 
if the system had been observed near periastron, it would have 
been difficult to confirm its binary nature. 



Comparison systems: Some binary systems have been identified 
where both components are more massive than those of R139: 
NGC-3603-A1, is a system comprised of a 1 16M0 primary and 
a 89M star ( |Schnurr et al.|2008a|>, and also WR20a, an 83M0 
and 82M system ( |Rauw et al.||2004| |Bonanos et al.| |2004). 
However, there are not many systems with a pair of massive 
evolved O-stars. Closer analogs are the Cyg OB2-B17 system 
( Stroud et al.|2010" l, where the component stars are 07 and 09 
supergiants and Cyg OB2-#5 where one of the stars is an 06-7 
supergiant (Ra uw et al.||1999) . It would appear that neither of 
these systems contain stars as massive as those predicted here. 
Consequently, it can be argued that R139 is the most massive O 
supergiant binary system yet discovered. 

Evolutionary masses: Figure [5] shows how the R139 system 
compares to evolutionary tracks from Friedrich et al. (in prep.), 
computed analogously to the models of Brott et al. ( 2011) l. The 
figure shows that the total luminosity and T e g of the system 
equals that of a single star with an initial mass above 125 Mq. 



Assuming R139 consists of coeval stars with a mass ratio 
of 1.2, the estimates for the current masses are 75 + 14Mq for 
the (cooler) primary and 62 + IIM0 for the (hotter) secondary. 
These values have been derived using a^- 2 method to fit the com- 
bined luminosity of the stars against that quoted for the system. 
The effective temperatures of the stars were fitted against the 
CMFGEN-derived temperatures of 34+2 kK. Interestingly, these 
estimated masses are in close agreement with the lower-mass 
limits from the orbital solution. This implies that the system has 
a high inclination and supports the need for additional photomet- 
ric observations. 

In these models an initial equatorial velocity of 110 kms -1 
was adopted in agreement with the current observed v sin i. The 
effect of rotation on the evolutionary tracks is very limited fo r 
initial rotation rates up to about 300 km s -1 ( Brott et al.||201 



Nevertheless, these tracks are sensitive to uncertain physical pro- 
cesses such as internal mixing and mass loss. The errors on the 
mass estimates represent the formal 1-sigma confidence limits 
of the x 2 fit and do not include systematic uncertainties in the 
model physics. 

The best fit corresponds to an age of 2 - 2.5 Myr and implies 
that both stars have significantly evolved off the zero-age main 
sequence. As the stars are assumed to be coeval, the substan- 
tial mass ratio implies a large difference in temperature between 
the component s (see Fig. [5]). T his is surprising given the similar 
spectral types; Martins et al. ( 2005 ) predict a temperature dif- 
ference nearer to 1 kK for a 0.5 variation in spectral types. This 
discrepancy may well reflect our still limited understanding of 
the physics of the most massive stars, illustrating the potential of 
massive binaries as tools to evaluate our models. 

The high quality, time-sampled VFTS observations have re- 
vealed that R139 is a binary system. The data suggest that it is 
the most massive evolved O-star binary system yet discovered: a 
result which further observations around periastron would help 
to confirm. As demonstrated here, such a massive system has al- 
ready presented challenges for theoretical models to reproduce 
its observed properties and it will likely provide a crucial test for 
evolutionary and atmospheric models in the future. 
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